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The conductance of a junction between a normal metal and a superconductor having the symmetry pro- 
posed by Berezinskii is studied theoretically. The main feature of this symmetry is the odd frequency depen- 
dence of the anomalous Green function, which makes possible the s-wave triplet superconducting state (the 
Berezinskii superconductor). The Andreev reflection (which links positive and negative energies) is sensitive to 
the energetic symmetry; as a result, the conductance of the junction involving the Berezinskii superconductor 
is qualitatively different from the case of a conventional superconductor. Experimentally, the obtained results 
can be employed to test the possibility of the Berezinskii superconductivity proposed for Na^CoO^ and to 
identify the odd-tJ component predicted for superconductor-ferromagnet systems. 
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The superconducting pairing can be described by 
the anomalous Green function F(l;2) (also referred 
to as the pair wave function) in the Matsubara tech- 
nique. The Pauli principle requires antisymmetry un- 
der permutation of two electrons in a Cooper pair, 
F(2; 1) = —F(l; 2), which leads to the standard classi- 
fication of superconducting phases [T] : if the coordinate 
dependence of F is even (s-wave, d-wave, etc.), then the 
spin dependence must be odd (singlet) and vice versa. 
This assumes that F is an even function of the imag- 
inary time t = 7~i — T2 or of the frequency u> in the 
Fourier representation. 

In 1974, Berezinskii suggested [2] that the s-wave 
triplet phase (not listed above) is also possible if F 
is an odd function of uj. The Pauli principle is satis- 
fied and there are no symmetry restrictions on the exis- 
tence of such a phase. Of course, the question remains 
whether an interaction necessary for the realization of 
the odd-w phase exists in real materials. Berezinskii 
discussed the odd-w phase as a possible phase of 3 He 
and argued (without a full microscopic derivation) that 
it could be formed due to the retarded paramagnon ex- 
change. However, up to now there are no indications for 
such a state in 3 He. Moreover, there is no microscopic 
theory producing the Berezinskii state in a bulk mate- 
rial. At the same time, properties of hypothetic odd-w 
superconductors were studied in a number of papers, 
e.g. ®. 
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Surprisingly, in 2001 the "long-range" Berezinskii 
superconductivity was theoretically discovered by Berg- 
eret, Volkov, and Efetov [4j (see also [3]) in a SF sys- 
tem consisting of a conventional (s-wave singlet) super- 
conductor and a ferromagnet. They demonstrated that 
the superconducting component with the symmetry pro- 
posed by Berezinskii arises in the case of inhomogeneous 
magnetization of the ferromagnet due to the proxim- 
ity effect and penetrates the ferromagnet over a much 
longer distance compared to the singlet component, 
since the exchange field does not suppress the triplet 
superconducting correlations with projections S z = ±1. 
This feature enables to spatially single out the Berezin- 
skii component. Many unusual properties of this be- 
havior have been recently investigated [6]. At the same 
time, experimental verification of the S z — ±1 compo- 
nent in SF structures is still under debate. A number 
of experiments observed a long-range proximity effect 
[7] and a Josephson coupling through a half-metallic 
ferromagnet [8], which can be interpreted in terms of 
the long-range component. However, the data on the 
structure of the magnetic inhomogeneity in experimen- 
tal samples is lacking. 

At the same time, as argued in |5J, the Berezinskii su- 
perconductivity can be realized in Na^OoO^. This con- 
jecture is based on the band-structure calculations and 
on the experimental results which indicate the triplet 
superconductivity (from the Knight-shift data) and the 
s-wave symmetry (from insensitivity to impurities). 

The question arises: can we suggest an experiment 
which is sensitive to the most nontrivial feature of the 
Berezinskii superconductivity, the odd-w dependence? 
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If yes, then this experiment could be a good test for 
such a state, similarly to the famous experiments sen- 
sitive to the nontrivial spatial symmetry of anisotropic 
superconductivity (the idea proposed in [lOj was em- 
ployed to experimentally verify the <i-wave symmetry 
in YBaCuO, see review [TT]). It seems natural that the 
odd-Lu superconductivity should lead to unusual features 
of the Andreev reflection [12] from such a superconduc- 
tor, since this process links an electron with positive 
energy to a hole with negative energy (with respect to 
the Fermi level). 

In this Letter, the differential conductance of the 
normal-metal — superconductor (NS) junction shown in 
Fig. [T] is studied at zero temperature. I consider three 
possibilities for the superconducting reservoir: S — a 
conventional superconductor with a gap (in this case 
previous results are reproduced, see reviews [T3] and 
references therein), So — a conventional superconduc- 
tor without a gap (e.g., due to paramagnetic impurities 
|14j). and Sb — a Berezinskii superconductor. Sb has 
two main features: odd-u symmetry and gapless spec- 
trum. The case of So is considered in order to reveal 
the features of Sb which are due to the unusual symme- 
try and not simply due to its gapless spectrum. Phys- 
ical examples of the Sb state are the bulk Berezinskii 
superconductor and the SF system (in the latter case, 
the normal wire should be attached to a region where 
only the long-range triplet component survives, while 
the short-range ones are negligible). 

The differential conductance of the junction is 



Gns{V) 



its normal-state value is Go 



dV 
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The 



diffusive normal wire of length L is characterized by the 
Thouless energy Et\i — D/L 2 , where D is the diffu- 
sion constant. At low voltages, the Andreev reflection 
plays an important role in the transport. I assume eV, 
i?Th "C Eo, where Eq is the energy scale on which the 
Green function of the superconductor varies (in a con- 
ventional superconductor Eq coincides with the static 
order parameter A). At the same time, the relation 
between eV and £/rh is arbitrary. 

I consider the dirty limit and employ the Usadel 
equation [15l [16] for the Green function g which is a 
8x8 matrix in the Keldysh ® Nambu-Gor'kov ® spin 
space: 

-PS- 

The Usadel equation for g is 



DV(gVg)+[iET 3 a i,g] =0, 



(3) 
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FigfT] NS junction. The conductances of the normal 
wire and the NS interface are Gn and Gt, respectively. 
The voltage V is applied to the normal reservoir N rea . 

where f and o denote the Pauli matrices in the Nambu- 
Gor'kov and spin spaces, respectively, while 1 is the unit 
matrix in the Keldysh space. At the interface with the 
normal reservoir (x — —L), g must be equal to the 
bulk normal-metallic function, while the Kupriyanov- 
Lukichcv boundary condition [17] at the NS interface 
(x = 0) reads 



Gn 



2GnL 



[g, as] 



(4) 



where gs is the Green function in the superconductor. 
The current is I = J jdE with the spectral current 



j(E, x) = Tr [T 3 a (gVg) 

loe 



A'l 



(5) 



Due to the normalization condition g 2 = 1, we can 
parametrize g K as g K = g R f — fg A , then the general 
relation g A = —t 3 (<7 fl ) ^3 allows us to consider only 
g R and / as independent functions. 

In S and So, we can write 



g R = T 3 a G R + n&oF 1 



(6) 



while Sb is a spin-triplet superconductor and the 
anomalous part is a vector in the spin space. Choos- 
ing its direction as z, we write [18] 



9 R 



f 3 a G R + f^F 1 



(7) 



Physically, this is a triplet superconducting state with 
zero projection of the Cooper pair's spin on the z axis, 
while the 1 and —1 projections on any axis in the xy 
plane are equiprobable. 

The Berezinskii superconductivity is odd in the 
Matsubara frequency, F(—qj) = —F(ui), therefore in 
the real-energy representation we obtain F R (—E) = 
—F A (E). Together with the general relation F A (E) = 
(F R (E)Y this yields Qjj] F R {-E) = - (F R (E))* . We 
need to know the low-energy behavior of F R (E), which 
depends on the choice of a model. The question of a mi- 
croscopic model for a bulk Berezinskii superconductor 
is not clear at present. At the same time, we know that 
such a state is realized in inhomogeneous SF structures 
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due to the proximity effect. Taking the low-energy be- 
havior of the odd-uj triplet component from, e.g., [20] . 
we obtain F(u>) = iAsgnui, hence F R (E) — iA at E — > 
0, with real constant A. This behavior also follows from 
the standard relation F R = A/VA 2 — E 2 in the case of 
the linear low-energy behavior A(E) = E/(l + A~ 2 ). 

The constraint (g R ) 2 = 1 allows us to parametrize 
the normal and anomalous Green functions by a (com- 
plex) angle as G R = cos and F R — sin 9 in the cases 
of both the conventional and Berezinskii superconduc- 
tor, Eqs. ([6]) and J7J. The equation and the boundary 
conditions for 9 are 



— 6" + iE sinO = 0, 



= 



x=—L 



Gt 
GnL 



sin 



-8) 



x=0 



(8) 
(9) 



The superconductor is described by 9s . In order to con- 
sider low voltages, we need to find the low-energy solu- 
tion of the Usadel equation. If E <C Eq, then s = ir/2 
in S, < 9 S < it/2 in So, and 9 S = id s with real d s 
in Sb- Thus the type of the superconductor enters our 
consideration only through the low-energy value 9s- 

In the case of conventional superconductors, S or So, 
components of the distribution function can be chosen 
as 

/ = tooo/o + T 3 & f 3 , (10) 

while in the case of Sb the number of components is 
doubled in the general case: 



/ = TQ&ofo + T 3 a f 3 + f <7 3 /o + T 3 (T 3 f 3 



(11) 



The spectral current ([5]) at x = can be rewritten 
with the help of the boundary condition ^ as 



j(E) = (G T /8e)/ 3 x 

< [{G R - G A )(G§ - Gg) + {F R + F A ){Fg + F#)] 



(12) 



The two terms in the square brackets have clear phys- 
ical meaning. Since (G R — G A )/2 is the single-particle 
density of states (DOS), the first term in the square 
brackets describes the quasiparticle contribution to the 
current. At the same time, the second term (of FF 
type) describes the supercurrent due to the Andreev re- 
flection. 

The subsequent derivation is similar to the conven- 
tional case [2lJ 03] • At zero temperature the integral 
over energies, which determines Gns{V), reduces to the 
sum of the two terms with E = ±eV. Finally, using the 
symmetry 9(—E) — 9*(E) for the conventional even-w 



superconductor or 
[15] , we obtain 

1 1 



-E) = -9*{E) for the odd-w one 



Gns(V) GnL J_ l cosh 9 2 (x) 



dx 



1 



Gtcos(9si — #i) cosh 0,52 cosh 02 



, (13) 



where the right-hand side (r.h.s.) is taken at E = eV 
and we have separated into the real and imaginary 
parts, = 0i + 202. The angles and 0s in the second 
term are taken at the NS interface. The problem now 
reduces to solving Eqs. ((8j- and calculating the r.h.s. 

of Eq. eng. 

We start from the simplest case of zero bias V = 0. 
At E = 0, Eqs. ||HJ) and ([9|) are solved by a linear func- 
tion. In the case of S or So we obtain 



G 



Grp 



N 



COS(#s 



(14) 



where 0o must be determined from the equation O 
(Gt/Gn) sin(0s — 0o)- In the case of Sb: 



n-i (n\ - G Jv ltanh ^o 



Grp 



0o 



cosh $o coshes ' 



(15) 



where $q must be determined from the equation i?q = 

An immediate consequence of these results is that 
G/vs(0) < Go in the cases of S and So, while Gats(O) > 
Go for Sb- In the S case, where the low-energy DOS is 
zero, the current is entirely due to the Andreev contri- 
bution. In the So case, a finite DOS appears and the 
current has both the Andreev and quasiparticle contri- 
butions. Interestingly, in the Sb case, only the quasi- 
particle processes contribute to Gats(O); this fact can 
be interpreted as the absence of the Andreev reflection 
from Sb at E — > 0. 

At eV <C i^Th i we can develop a perturbation theory 
finding corrections to the zero-bias solution of the Us- 
adel equation. Two orders in eV/-&rh yield a quadratic 
low-bias behavior: Gns(V) = Gjvs(O) + aV 2 . The ex- 
plicit form of a is cumbersome and I only present its 
most important features. The sign of a depends on 
the ratio Gn/Gt and on the type of the superconduc- 
tor: a > for S and So if Gn/Gt < g c and a < 
if Gn/Gt > g c where g c is of the order of unity and 
weakly depends on 05 (as 9s changes from tt/2 to 0, the 
critical value g c stays between 0.8 and 0.9), while a < 
for Sb at arbitrary Gn/Gt- 

Now let us consider in more detail the limiting cases 
of large and small ratio Gn/Gt- 
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Tunneling limit: Gn ^ Gt. In this case, we can 
retain only the second term in the r.h.s. of Eq. (fT3")) . 
The proximity effect is weak, \9\ <C 1. 

In the case of Sq and Sb we can set 9 = 0, then 



G NS (V) 
Go 



(16) 



where vs{E) is the DOS in the superconductor, normal- 
ized to the normal-metallic value; this is nearly constant 
at E <C Eg (the more accurate analysis presented above 
gives the small bias-dependent correction to this con- 
stant). The physical meaning of Eq. ([To]) is the tunnel- 
ing spectroscopy of the superconductor with the normal 
probe. An essential difference between Sq and Sb is that 
^s(O) = cos 9s < 1 for So, while ^s(O) = coshes > 1 
for Sb (this fact for Sb was pointed out in [T5]). 

In the case of S, us(eV) = below the gap, there- 
fore we cannot neglect the proximity effect. Linearizing 
Eqs. and © over 9, we find the solution and finally 
obtain 

G NS {V) _G T sinh(2 A /i) + sin(2Ve) eV 

Go ~ G N 4^[sinh 2 ( y/e) + cos 2 {y/e)\ ' 6 ~ iW 

(17) 

hence Gns(V) <C Go [Eq. (fTF)) also follows from a more 
general result obtained in 21j]. 

The results for the tunneling limit arc summarized 
in Fig. [2j Note that the physics related to the Andreev 
reflection is not essential for Sb in this limit, since the 
transport is due to the quasiparticle contribution. At 
the same time, the Andreev reflection plays a crucial 
role in the transparent limit considered below. 

Transparent limit: Gn <C Gt- In this case, we 
can retain only the first term in the r.h.s. of Eq. (|13p . 

At eV <C -Exh we calculate a small correction to the 
zero-bias conductance. For S and So we obtain 



G NS {V) 
Go 



= 1 + A 



( eV 



\-Erh 



with the positive coefficient 



A 



sin 2 9 S 3 sin 29 s 2 sin 2 9 S 



49 s 



(18) 



(19) 



A(9s) monotonically grows from zero at 9s = to ap- 
proximately 0.015 at 9s = 7r/2. 

The low-bias conductance for Sb is 



Gns(V) 



c 



Go tanh s 

with the positive coefficient 
9 3 

B = 



B 



eV 
~E' 



Th 



(20) 



3cosh 2 S 4sinh20 s 



0| tanh 2 0c 



0| tanh 5 
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(21) 
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ST 
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eV/E^ 



Fig[2] Differential conductance at eV <C Eq in 
the tunneling limit (Gn S> Gt). In the S case, 
Gns(V) demonstrates the zero- bias anomaly [22] : it 
has quadratic and 1/VV dependences at eV <C i?Th 
and eV S> -&rh, respectively. In the So and Sb cases, 
Gns(V) is nearly constant, smaller than Go for So and 
larger than Go for Sb- 




eV/E Th 

FigE] Differential conductance at eV <C -Eo in the 
transparent limit (Gn <C Gt). In all the three cases, 
Gns(V) is quadratic at eV <C ETh and approaches 
unity as 1/y/V at eV 3> ETh- In the S and So cases, 
the behavior is reentrant (see [23] for the S case). The 
maximum value of Gns(V)/Go (achieved at eV of the 
order of several Exh) approximately equals 1.15 for the 
S case, while for So the curve is closer to unity. In the 
Sb case, Gjvs(V) monotonically decreases. 



B^s) is a monotonically growing function starting from 
zero at 0s =0. Although 05 is an unknown parameter, 
it can in principle be determined from measurements in 
the tunneling limit [see Eq. (fT6"j) with z^s(0) = cosh0g]. 

At eV ^> E^ all the three cases (S, So, and Sb) 
are treated in a similar manner. Since Gn <C Gt, the 
boundary condition JSJ at x = reduces to 9 = 9s- 
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The well-known solution of the sine-Gordon equation 
dS]) with fixed surface value is 

0(x) — 4arctan|tan {^~^J ex P^ — (1 — ^ a 'l\/^) V 

(22) 

This solution satisfies the boundary condition 6{—L) = 
with good accuracy, because 8(—L) is exponentially 
small at E E Th . Substituting Eq. (12"2")) into the first 
term in the r.h.s. of Eq. (fT5|) , we obtain 

(23) 

where the positive coefficient C depends only on 9s, 
i.e., on the type of the superconductor. In the S case, 
C«0.3. 

The results for the transparent limit are summa- 
rized in Fig. [3] Note that in 0] and [B] a different, 
so-called cross geometry was considered under assump- 
tion of weak proximity effect. Then a small correction 
to the normal-state conductance due to the Sb compo- 
nent was numerically shown to monotonically decrease 
as a function of temperature at zero bias. 

In conclusion, the conductance of the junction be- 
tween a normal metal and a Berezinskii superconduc- 
tor (odd-w spin-triplet s-wave state) has been studied. 
The main differences from the case of a conventional 
superconductor are: (i) in the tunneling limit, Gns(Y) 
is larger than the normal-state conductance (Fig. [2]), 
and (ii) in the transparent limit, Gns(V) monotonically 
decreases (Fig. [3]). These features can be used as an 
experimental test for the Berezinskii superconductivity 
in bulk samples (e.g., Na^OoO^) or in superconductor- 
ferromagnet proximity systems. 
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